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Abstract 

We have re-analyzed the world data on inclusive polarized DIS, in both NLO 
and LO QCD, including the new HERMES and COMPASS data. The updated 
NLO polarized densities are given in both the MS and JET schemes. The impact 
of the new data on the results is discussed. 
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Since the famous European Muon Collaboration (EMC) experiment [1] at CERN 
in 1987, substantial efforts, both experimental and theoretical, have been devoted to 
understanding the partonic spin structure of the nucleon, he., how the nucleon spin is 
built up out from the intrinsic spin and orbital angular momentum of its constituents, 
quarks and gluons. Our present knowledge about the spin structure of the nucleon 
comes mainly from polarized inclusive and semi-inclusive D1S experiments at SLAC, 
CERN, DESY and JLab, polarized proton-proton collisions at RHIC and polarized 
photoproduction experiments. The determination of the longitudinal polarized parton 
densities in QCD is one of the important and best studied aspects of this knowledge. 

In this Brief Report we present an updated version of our Set 1 and Set 2 NLO 
QCD polarized parton densities in both the MS and the JET (or so-called chirally 
invariant) [2] factorization schemes, as well as the LO ones, determined in our recent 
analysis [3]. Comparing to our previous analysis: i) The old HERMES proton [4] and 
neutron [5] data are replaced with the new HERMES proton and deuteron data [6] and 
ii) a complete treatment of the recent COMPASS data on the longitudinal asymmetry 
Af [7] is included. A FORTRAN package of the obtained polarized parton densities 
(LSS’05) will be presented at the Durham HEPDATA web site to be used for practical 
purposes. 

In QCD the spin structure function r/j has the following form ( Q 2 >> A 2 ): 

9i( x , Q 2 ) — 9i( x -> Q 2 ) lt + 9i(x, Q 2 ) HT , (1) 

where ”LT ,? denotes the leading twist (r = 2) contribution to g\ , while ”HT” denotes 
the contribution to g\ arising from QCD operators of higher twist, namely r > 3. In 
Eq. (1) (the nucleon target label N is dropped) 

9i( x , Q 2 ) LT = <71 (X, Q 2 ) p qcd + h™ c (x, Q 2 )/Q 2 + 0(M 4 /Q 4 ) , (2) 

where gi(x, Q 2 ) p qcd is the well known (logarithmic in Q 2 ) pQCD contribution and 
h TMC (x,Q 2 ) are the calculable [8] kinematic target mass corrections, which effectively 
belong to the LT term. In Eq. (1) 

9 i(x, Q 2 ) HT = h(x, Q 2 )/Q 2 + 0(A 4 /Q 4 ) , (3) 

where h(x, Q 2 ) are the dynamical higher twist (r = 3 and r = 4) corrections to <?i, which 
are related to multi-parton correlations in the nucleon. The latter are non-perturbative 
effects and cannot be calculated without using models. 
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Let us recall that the Set 1 polarized parton densities correspond to fits to gi/Fi 
and Ai(pz g 1 /F\) data (so called ’ gi/Fi fits): 


gi{x,Q 2 ) A ^i(^,Q 2 ) p qcd 

Fi(x,Q*)\ exp Fi(*,g 2 )pQCD ’ 


( 4 ) 


where for the structure functions g\ and F\ their leading twist NLO QCD expressions 
are used. (Why this method is incorrect for extracting the LO polarized PDs, is 
discussed in Ref. [11].) The Set 2 polarized PD correspond to fits to g\/F\ and A 1 
data where the experimental data for the unpolarized structure function Fi(x,Q 2 ) are 
used „ „ 

g?(x,Q 2 ) LT + h N (x)/Q 2 


g?(x,Q 2 ) 

Fi*(x, Q 2 ) 




J exp 


Fi {x, Q 2 ) 


( 5 ) 


exp 


As usual, Fi is replaced by its expression in terms of the usually extracted from 
unpolarized DIS experiments F 2 and R and phenomenological parametrizations of the 
experimental data for F 2 (x, Q 2 ) [9] and the ratio R(x, Q 2 ) of the longitudinal to trans¬ 
verse 'yN cross-sections [10] are used. Note that such a procedure is equivalent to a fit 
to gi data themselves and we will refer to these as \g\ T + HT)’ fits. In this case the 
HT corrections to g\ cannot be compensated and have to be taken into account (for 
more details see our previous paper [3]). In (5) gi(x, Q 2 )lt (N=p, n, d) is given by 
the leading twist expression (2) in LO/NLO approximation including the target mass 
corrections and h N (x) are the dynamical r = 3 and r = 4 HT corrections which are 
extracted in a model independent way. 

As in our previous analyses [3, 12], for the input LO and NLO polarized parton 
densities at Q‘q — 1 GeV 2 we have adopted a simple parametrization 

xAu v (x,Q 2 0 ) = g u A u x au xu v (x,Ql), 
xAd v (x,Ql) = g d A d x ad xd v (x,Ql), 
xAs(x,Ql) = g s A s x as xs(x,Q 2 0 ), 

xAG(x, Ql) = g g A g x a9 xG(x,Ql), (6) 

where on the RHS of (6) we have used the MRST98 (central gluon) [13] and MRST99 
(central gluon) [14] parametrizations for the LO and NLO (MS) unpolarized densities, 
respectively. The normalization factors Ai in (6) are fixed such that rj l are the first 
moments of the polarized densities. To fit better the data in LO QCD, an additional 
factor (1+7 j,x) on the RHS is used for the valence quarks. Bearing in mind that the light 
quark sea densities A u and Ad cannot, in principle, be determined from the present 
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inclusive data (in the absence of polararized charged current neutrino experiments) we 
have adopted the convention of a flavor symmetric sea 

A u sea = A u = A d sea = Ad = As = As. (7) 

Note that this convention affects the results for the valence parton densities, but not 
the results for the strange sea quark a gluon densities. 

The first moments of the valence quark densities r/„ and r]d are constrained by the 
sum rules 

a 3 = gA = F + D = 1.2670 ± 0.0035 [15], (8) 

a 8 = 3F-D = 0.585 ± 0.025 [16], (9) 

where a 3 and a 8 are non-singlet combinations of the first moments of the polarized 
parton densities corresponding to 3 rd and 8 th components of the axial vector Cabibbo 
current 

a 3 = (A u + A u)(Q 2 ) - (Ad + Ad)(Q 2 ) , (10) 

a 8 = (A u + Au)(Q 2 ) + (Ad + Ad)(Q 2 ) - 2(As + A s)(Q 2 ) . (11) 

As in [3], we have used the MRST02(NLO) unpolarized parton densities [17] to 
constrain via positivity our Set 1 and Set 2 polarized PD in both MS and JET schemes. 

The numerical results of our fits to the world data [1, 6, 7, 18] on cj\/F\ and A x 
are presented in Tables I and II. The data used (190 experimental points) cover the 
following kinematic region {0.005 < x < 0.75, 1 < Q 2 < 58 GeV 2 }. The total 

(statistical and systematic) errors are taken into account. The systematic errors are 
added quadratically. It is seen from the Tables I and II that the values of y 2 /DF(MS) 
and y 2 /DF(JET) coincide almost exactly for the 'g\/F\ as well as for the (g\ T + HT) 
fits, which is a good indication of the stability of the analysis regardless of the scheme 
used. Analytic formulae for the input parton densities and higher twist contributions 
are given in the Appendix. 

Let us note the main features of our new results: 
i) The new Set 1 and Set 2 polarized PD are close to those determined in our 
recent analysis [3]. However, due to more accurate HERMES/d and COMPASS (small 
x region) data the polarized PD are better determined now (see the errors in the Tables 
1 and 2 in [3] and this paper). 

ii) Compared to our previous results [11, 12] we obtain now smaller values for the 
gluon polarization (the first moment of A G(x, Q 2 )), which leads to a smaller difference 
between the values of the strange quark polarization (the first moment of As(x, Q 2 )) 
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determined in the MS and the JET schemes, respectively. As a consequence, the 
difference between the quark helicity AS determined in these two schemes is also 
smaller. This fact has been already observed in our recent analysis [3]. 

Hi) The shape of both the polarized strange quark and gluon densities is rather 
different from that of LSS’Ol [12] (see Fig. 1). This is mainly due to the different 
positivity bounds imposed on the polarized PD in the LSS’Ol and LSS’05 fits. The 
impact of positivity constrains on polarized PD has been discussed in detail in [3]. 




Figure 1: Comparison between our two sets of NLO(MS) polarized parton densities, 
LSS’Ol and LSS’05(Set 1), at Q 2 = 4 GeV 2 . 

iv) The effect of the new data on the higher twist corrections to the proton and 
neutron spin structure functions, h p (x ) and h n (x), is negligible (see Fig. 2). The 
new values are in good agreement with the old ones although there is a tendency for 
the central values for the proton target to be slightly higher than the old ones. For 
the neutron target the only difference is that the new value of h n (x) at x ~ 0.2 is 
considerably higher than the old one and definitely different from zero. 

Finally, let us illustrate how important the higher twist corrections are for the 
description of the data. In Fig. 3 we compare the new very accurate HERMES deuteron 
data on the ratio ()\/F\ at measured x and Q 2 with the theoretical curves for g x obtained 
by: i) the best fit to the data using only the LT term in Eq. 5 (dotted curve) and ii) 
the best fit to the data taking into account in (5) the HT corrections too (solid curve). 
The dashed curve corresponds to the LT term when the HT corrections are accounted 
for. As expected, the role of HT corrections is significant in the small x region where 
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the values of Q 2 are relatively small: Q 2 rs 1.2 — 2.5 Gel/ 2 . 


NLO MS 



Figure 2: Effect of the new data on the higher twist values. 



x 


Figure 3: The effect of higher twist corrections (see the text). 
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Appendix 


A. Input Parton densities 

For practical purposes we present here explicitly our Set 1 and Set 2 of polarized 
parton densities at Q 2 = 1 GeV 2 . The polarized valence quark densities correspond to 
the convention of a SU(3) flavour symmetric sea. 


LSS’05 (Set 1) - NLO(MS) PD^/Fi): 


xA u v (x) 

= 0.4604 x a6240 (1 - a;) 3 ' 428 ( 

1 + 2.179 x 1/2 + 14.57 x ) , 


xA d v {x) 

= -0.02408 a; 0 ' 3680 (1 - a;) 3 ' 864 

( 1 + 35.47 x l/2 + 28.97 x ) , 


xAs(x) 

- -0.02226 a; 0 ' 3546 (1 - a;) 7 ' 649 

( 1 + 3.656 x 1/2 + 19.50 x ) , 


xAG(x) 

= 641.0 a; 3 ' 341 (1 - a:) 6 ' 879 ( 1 

- 3.147 x 1/2 + 3.148 x ) . 

(A.l) 

LSS’05 (Set 1) - NLO(JET) PD(^/Fi): 



xA u v {x) 

= 0.4589 a; 0 ' 6224 (1 - x) 3A28 ( 

1 + 2.179 x 1/2 + 14.57 x ) , 


xA d v {x) 

= -0.02490 a; 0 ' 3813 (1 - a;) 3 ' 864 

( 1 + 35.47 x 1/2 + 28.97 x ) , 


xAs(x) 

- -0.01577 x 0 - 31 ' 27 (1 - a;) 7 ' 649 

( 1 + 3.656 x 1/2 + 19.50 x ) , 


xAG(x) 

= 923.6 a; 4 ' 138 (1 - a;) 6 ' 879 ( 1 

- 3.147 x 1/2 + 3.148 x ) . 

(A.2) 

LSS’05 (Set 2) - LO PD(^i + HT): 



xAu v (x) = 

0.1760 a; 0 ' 3012 (1 - x) 3A77 (1 + 1.610 a;)( 1 - 0.4085 x 1/2 + 17.60 

x ) , 

xA d v (x) = - 

-0.00765 a; 0 ' 1535 (1 - a:) 3 ' 398 (1 + 3.797 a;)( 1 + 37.25 x l/2 + 31.14 

x ) , 

xA s(x) = - 

-0.04660 a; 0 ' 3542 (1 - a;) 8 ' 653 ( 1 - 0.9052 x l/2 + 11.53 x ) , 


xA G(x) = 

303.6 a: 3 ' 188 (1 - a;) 6 ' 879 ( 1 - 3.147 x l/2 + 3.148 x ) . 

(A.3) 

LSS’05 (Set 2) - NLO(MS) PD(cq + HT): 



xA u v (x) 

= 0.5041 a; 0 ' 6689 (1 - x) 3A 28 ( 

1 + 2.179 x 1/2 + 14.57 x ) , 


xA d v (x) 

= -0.02847 x 0A364 (1 - a;) 3 ' 864 

( 1 + 35.47 x 1/2 + 28.97 x ) , 


xAs(x) 

= -0.02759 x 0 - 3740 (1 - x) 7 ' 649 

( 1 + 3.656 x 1/2 + 19.50 x ) , 


xAG(x) 

= 303.6 a: 3 ' 188 (1 - a:) 6 ' 879 ( 1 

- 3.147 x 1/2 + 3.148 x ) . 

(A.4) 
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LSS’05 (Set 2) - NLO(JET) PD(^i + HT): 

xAu v (x) = 0.4991 x a6639 (1 - a:) 3 ' 428 ( 1 + 2.179 x l/ 2 + 14.57 x ) , 

xAd v (x) = -0.02693 a; 0 ' 4133 (1 - a;) 3 ' 864 ( 1 + 35.47 x 1/2 + 28.97 x ) , 
xAs(x) = -0.02434 a; 0 ' 3837 (1 - a;) 7 ' 649 ( 1 + 3.656 x l/2 + 19.50 x ) , 
xAG(x) = 736.1 a: 3 ' 678 (1 - a;) 6 ' 879 ( 1 - 3.147 x l/2 + 3.148 x ) . (A.5) 

B. Higher twist corrections 

We present here the best fit parametrizations of the extracted values of higher twist 
corrections to g\ for the proton and neutron targets (see Table 2) which are valid in 
the experimental x region: 0.005 < x < 0.75 . 

Fit \g±(x, <5 2 ){lo+tmc} + h(x)/Q 2 ]: 


h p (x ) 

0.2186 r „ . 

- 0.0223 __ exp[ 2({x 0.1905)/0.1761) 2 ] 



V 77 / 2 


h n (x ) 

0.2428 r .. w 

- 0.0362+ _ exp[ 2{(x 0.0267)/0.1768) 2 ] 

(B.l) 


V 7r /2 


Fit \gi(x, <5"){nlo(ms)+tmc} + h{x)/Q 2 \. 

n i '■70 (-\ 

h p (x) = 0.0360 - '1— exp[-2((x - 0.2030)/0.1736) 2 ] 

0 2313 

h n (x ) = 0.0176 H— '-j^=exp[—2((x — 0.0723)/0.1748) 2 ] (B.2) 

H 2 

Fit \gi(x, Q 2 ){nlo(jet)+tmc} + h{x)/Q 2 \. 

0 1629 

h v (x) = 0.0373 - ‘_ ea^[-2((a; - 0.2046)/0.1700) 2 ] 

H 2 

0 2519 

h n (x) = 0.0137 H— '-j^=exp[—2{{x — 0.0676)/0.1812) 2 ] (B.3) 

H 2 
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TABLE I. The parameters of the Set 1 NLO input polarized PD at Q 2 = 1 GeV 2 
as obtained from the best 'gl/FV fits to the world data in the MS and JET schemes. 
The errors shown are total (statistical and systematic). The parameters marked 
by (*) are fixed by the sum rules (8) and (9). 


Fit 

r/ NLO/ F NLO (M g) 

^n lo / f n l °(jet) 

DF 

190 - 6 

190 - 6 

x 2 

166.2 

168.3 

X 2 /DF 

0.903 

0.915 

Vu 

0.926* 

0.926* 

&U 

0.207 ± 0.018 

0.206 ± 0.018 

Vd 

- 0.341* 

-0.341* 

&d 

0.098 ± 0.059 

0.111 ± 0.059 

Vs 

- 0.061 ± 0.007 

- 0.051 ± 0.008 

&S 

0.637 ± 0.082 

0.595 ± 0.101 

Vg 

0.307 ± 0.175 

0.171 ± 0.166 

dg 

2.372 ± 0.691 

3.169 ± 1.323 
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TABLE II. The parameters of the Set 2 of LO, NLO(MS) and NLO(JET) input 
polarized PD at Q 2 = 1 GeV 2 as obtained from the best (g\ + HT) fits to the 
world data. The errors shown are total (statistical and systematic). The 
parameters marked by (*) are fixed. Note that the TMC are included in {gi) hT . 


Fit 

(9i)Wo + h( X )/Q 2 

(^i)nIo(ms) + h (x)/Q 2 

(^i)nlo(jet) + h(x)/Q 2 

DF 

190 - 16 

190 - 16 

190 - 16 

x 2 

156.1 

154.5 

154.4 

x 2 /df 

0.897 

0.888 

0.887 

Vu 

0.926* 

0.926* 

0.926* 

CL U 

0.000 ± 0.010 

0.252 ± 0.037 

0.247 ± 0.038 

In 

1.610 ± 0.301 

0* 

0* 

Vd 

- 0.341* 

-0.341* 

-0.341* 


0.000 ± 0.062 

0.166 ± 0.124 

0.143 ± 0.149 

Id 

3.797 ± 1.736 

0* 

0* 

Vs 

- 0.068 ± 0.006 

- 0.070 ± 0.008 

- 0.060 ± 0.012 

CL S 

0.544 ± 0.058 

0.656 ± 0.069 

0.666 ± 0.114 

Vg 

0.179* 

0.179 ± 0.267 

0.231 ± 0.283 

dg 

2.218* 

2.218 ± 1.650 

2.709 ± 1.567 

Xi 

hP{xi) [GeV 2 ] 

0.028 

-0.010 ± 0.042 

0.018 ± 0.047 

0.021 ± 0.043 

0.100 

- 0.078 ± 0.034 

- 0.031 ± 0.032 

- 0.023 ± 0.034 

0.200 

- 0.151 ± 0.038 

- 0.100 ± 0.040 

- 0.093 ± 0.043 

0.350 

- 0.011 ± 0.044 

0.004 ± 0.046 

0.007 ± 0.046 

0.600 

0.022 ± 0.021 

0.036 ± 0.020 

0.037 ± 0.021 

Xi 

h n (xi) [GeV 2 ] 

0.028 

0.230 ± 0.061 

0.182 ± 0.065 

0.195 ± 0.063 

0.100 

0.174 ± 0.040 

0.196 ± 0.038 

0.203 ± 0.038 

0.200 

0.064 ± 0.054 

0.081 ± 0.061 

0.078 ± 0.065 

0.325 

0.038 ± 0.026 

0.025 ± 0.029 

0.023 ± 0.031 

0.500 

0.036 ± 0.014 

0.014 ± 0.013 

0.013 ± 0.015 
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